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LONG-TIME BEHAVIOUR OF A FULLY DISCRETE LAGRANGIAN 
SCHEME FOR A FAMILY OF FOURTH ORDER 

HORST OSBERGER 


Abstract. A fully discrete Lagrangian scheme for solving a family of fourth order equations 
numerically is presented. The discretization is based on the equation’s underlying gradient 
flow structure w.r.t. the L 2 -Wasserstein distance, and adapts numerous of its most important 
structural properties by construction, as conservation of mass and entropy-dissipation. 

In this paper, the long-time behaviour of our discretization is analyzed: We show that 
discrete solutions decay exponentially to equilibrium at the same rate as smooth solutions of 
the origin problem. Moreover, we give a proof of convergence of discrete entropy minimizers 
towards Barenblatt-profiles or Gaussians, respectively, using T-convergence. 


1. Introduction 

In this paper, we propose and study a fully discrete numerical scheme for a family of nonlinear 
fourth order equations of the type 

8 t u = — {u{u a ~ 1 u^ x ) x ) x + \{xu) x for x e D = R, t > 0 (1) 

and m(0, .) = u° on 12 at initial time t = 0. The initial density u° ^ 0 is assumed to be compactly 
supported and integrable with total mass M > 0, and we further require strict positivity of 
u° on supp(u°) = [a, b]. For the sake of simplicity, let us further assume that M = 1. We 
are especially interested in the long-time behaviour of discrete solutions and their rate of decay 
towards equilibrium. For the exponent in 0. we consider values a e [|, 1], and assume A ^ 0. 
The most famous examples for parabolic equations described by (jT|) are the so-called DLSS 
equation for a = (first analysed by Derrida, Lebowitz, Speer and Spohn in |24l 125] with 
application in semi-cunductor physics) and the thin-film equation for a = 1 indeed, for other 
values of a, references are very rare in the literature, except [45] of Matthes, McCann and Savare. 

Due to the physically motivated origin of equation (JT]) (especially for a = | and a = 1), it is 
not surprising that solutions to 0 carry many structural properties as for instance nonnegativity, 
the conservation of mass and the dissipation of (several) entropy functionals. In section[2j we are 
going to list more properties of solutions to ([!]). For the numerical approximation of solutions to 
0, it is hence natural to ask for structure-preserving discretizations that inherit at least some 
of those properties. A minimum criteria for such a scheme should be the preservation of non¬ 
negativity, which can already be a difficult task, if standard discretizations are used. So far, many 
(semi-)discretizations have been proposed in the literature, and most of them keep some basic 
structural properties of the equation’s underlying nature. Take for example mmmm, where 
positivity appears as a conclusion of Lyaponov functionals — a logarithmic/power entropy [10J 
1161 l40j or some variant of a (perturbed) information functional. But there is only a little number 
of examples, where structural properties of equation |l]) are adopted from the discretization by 
construction. A very first try in this direction was a fully Lagrangian discretization for the DLSS 
equation by During, Matthes and Pina [ 25 ], which is based on its L 2 -Wasserstein gradient flow 
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representation and thus preserves non-negativity and dissipation of the Fisher-information. A 
similar approach was then applied [47j . again for the special case a = where we even showed 
convergence of our numerical scheme, which was - as far as we know - the first convergence 
proof of a fully discrete numerical scheme for the DLSS equation, which additionally dissipates 
two Lyapunov functionals. 


1.1. Description of the numerical scheme. We are now going to present a scheme, which is 
practical, stable and easy to implement. In fact our dicretization seems to be so mundane that 
one would not assume any special properties therein, at the first glance. But we are going to show 
later in section [2j that our numerical approximation can be derived as a natural restriction of 
a L 2 -Wasserstein gradient flow in the potential landscape of the so-called perturbed information 
functional 


F a ,\{ v)=tt- f (d x u a ) 2 dx + ^ f \x\ 2 u{x)dx , 
2 a Ja 2 J n 


( 2 ) 


into a discrete Lagrangian setting, thus preserves a deep structure. The starting point for our 
discretization is the Lagrangian representation of 0- Since each u(t, •) is of mass M, there is 
a Lagrangian map X(t, •) : [0, M] —* Q — the so-called pseudo-inverse distribution function of 
u{t, ■) — such that 

r-x(U) 


e = 


1 


u(t, x) dx, for each £ e [0, M\. 


(3) 


Written in terms of X, the Wasserstein gradient flow for T a \ turns into an L 2 -gradient flow for 


1 f M 

T at x{uoX) = — 

2a J n 


X? 


-U 




rM 

1 : 


x 2 d£, 


that is, 
d t X = 


2o i i I 

( 2 a+l)* d *(Z a+ * d K Za+ *)+ AX, where Z(t,f) :=-^^ = u(t,X(t,f)). (4) 

To build a bridge from @ to the origin equation 0 , remember that 0 can be written as a 
transport equation, 

'SiF a ,x(u) 


dtu + (itv a ) x = 0, with velocity field v Q = — 


Su 


(5) 


where 5iF at \(u)/6u denotes the Eulerian first variation. So take the time derivative in equation 
(|3| and use ([5]), then a formal calculation yields 

rX(W 

0 = d t X(t, X(f, 0) + d t u(t, X(t, 0) dx 

J — 00 


J — OO 

r x(t,o 


rMh i) 

= d t X(t,€)u(t,X(t,€)) - (uv a )x(t,x) dx = d t X(t,£)u(t,X(t,€)) - (uv a )oX(^). 

J — 00 

This is equivalent to 

8 t X(t, 0 = v a o X(t, 0 for (t , 0 e (0, +oo) x [0, M], 


which is further equivalent to Q. 

Before we come to the proper definition of the numerical scheme, we fix a spatio-temporal dis¬ 
cretization parameter A = (t; 5): Given r > 0, introduce varying time step sizes t = (n, T 2 ,...) 
with r n e (0, r], then a time decomposition of [0, +oo) is defined by {t n } ns ^ with t n := Xq=i T n■ 
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As spatial discretization, fix K e N and S = M/K, and declare an equedistant decomposition of 
the mass space [0, M\ through the set {fk)k=o with := kS, k = 0,..., K. 

Our numerical scheme is now defined as a standard discretization of equation Q: 

Numerical scheme. Fix a discretization parameter A = (t; <5). Then for any (a, A) e [|, 1] x 
[0, +go) and any initial density function u° G L 1 (fl) satisfying the above requirements, a numer¬ 
ical scheme for 0 is recursively given as follows: 

(1) For n = 0, define an initial sequence of monotone values x^ := (a;[J, ... ,x° K ) e 
uniquely by xjj = a, Xk = b and 


= 


r 


u°( x) da;, for any k = 1,..., K — 1. 


The vector ~x° A describes a non-equidistant decomposition of the support [a, b] of the initial 
density function u°. In any interval [a;^_ 1; a;°], k = 1,..., K, the density u° has mass S. 
(2) For n > 1, define recursively a monotone vector x A := (xq, ..., x^) e M* r+1 as a solution 
of the system, consisting of (K + 1 )-many equations 

n— 1 


X, 


2a 


(^ +| r+§[D^+i] fc+ , +Xx k , (6) 


r„ (2a + 1) 2 <5 L v k+ ?‘ 

with k = 0,..., K, where the values zf_ i ^ 0 are defined by 

( _ s_ 

_ J x ] 

^ 0— 1 


0 


, for I = 1,..., K 
, else 


(7) 


and 


[D^] fc i := b~ 2 (z 


fc+i 


2~“+2 , “+2t 

Z "k- i +Z k-i)- 


We later show in Proposition \U\ that the solvability of the system 0 is guaranteed. 
The above procedure (1) — (2) yields a sequence of monotone vectors xa := (x^, x^, ..., x^,...), 
and any entry x^ defines a spatial decomposition of the compact interval [x^x^] c fl, n e N. 
Fixing k = 1,... ,K, the sequence n i-» x% defines a discrete temporal evolution of spatial grid 
points in Q, and if one assigns each interval [x^_ 1 ,x^] a constant mass package S, the map 
n i-» [x]^_ 1 ,x^] characterizes the temporal movement of mass. Hence xa is uniquely related to 
a sequence of local constant density functions ua '■= ..., u \,...), where each function 


: fi —» R+ holds 


K 


“a( x ) = U 5 [x^] := ^ 


=i J 'fc 


k = 1 


5 

x'u — X 


-I 


fc-1 


(x k - lt x k ] 


{x). 


( 8 ) 


We will see later in section 2.1 that the information functional F a ,\ can be derived using the 
dissipation of the entropy 

f \x\ 2 u(x) dx, with ip a (s) := 
z Jo (fcq^sliHs), a =2 

with constants 0 Q := i/2a/(2a + 1), and A a ,\ := •y / A/(2a + 1). As replacements for the entropy 
TLa t \ and the perturbed information functional F a ,\, we introduce 
k \ k ( 


H. a ,\{y) = <Pa(u) da; + 

Jn 


H n 


a(x) := <5 ^ fa Ofc-§) + 2 l^l 2 ’ with ) := 


k=1 


k =0 


a—1/2 > a ^ ( 2 i 1] 
0i/2 ln(s), OL=\ 


(9) 
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and 


F q> a(x) 


E 

kel° K 




<5El^l 2 - 


k =0 


( 10 ) 


1.2. Familiar schemes. The construction of numerical schemes as a solution of discrete Wasser- 
stein gradient flows with Lagrangian representation is not new in the literature. Many approaches 
in this spirit have been realised for second-order diffusion equation |9] HT] [44] |49], but also for 
chemotaxis systems [B], for non-local aggregation equations [1312D], and for variants of the 
Boltzmann equation j33 . We further refer [43] to the reader interested in a very general nu¬ 
merical treatement of Wasserstein gradient flows. In case of fourth order equations, there are 
some results for the thin-film equation and its more general version, the Hele-Shaw flow, see 
[5T1 331 . but converegence results are missing. Rigorous stability and convergence results for 
fully discrete schemes are rare and can just be found in |32J HB] for second order equations, and 
in m for the DLSS equation. However, there are results available for semi-discrete Lagrangian 
approximations, see e.g. H EZ| . 


1.3. Main results. In this section, fix a discretization A = (r; <5) with t,S > 0. For any 
solution xa of we will further denote by «a = («a,Ua) • • •) the corresponing sequence of 
local constant density functions, as defined in ([8]). 

All analytical results that will follow, arise from the very fundamental observation, that solu¬ 
tions to the scheme defined in section [ 13 ] can be successively derived as minimizers of the discrete 
minimizing movemend scheme 

+ ( n ) 

LTn k 

An immediate consequence of the minimization procedure is, that solutions x^ dissipate the 
functional F Qj a- 

Concerning the long-time behaviour of solutions xa , remarkable similarities to the continuous 
case appear. Assuming first the case A > 0, it turns out that the unique minimizer x^ 1111 of 
H q a is even a minimizer of the discrete information functional F a> A, and the corresponding set 
of density functions = u^Jx^ 1111 ] converges for 5 —» 0 towards a Barenblatt-profile b^A or 
Gaussian b]y 2 ,A, respectively, that is defined by 

b Q ,A = (a - &M 2 ) + (a 1/2) , b = i- /2 Aa,A if a. > 1/2 and (12) 

\ la 

bi/ 2 ,A = ae _Al / 2 ' x|:c|2 if a = 1/2, (13) 

where a e R is chosen to conserve unit mass. Beyond this, solutions Xa satisfying § converge as 
n —* oo towards a minimizer x)/ 1111 of F aj A with an exponential decay rate, which is ’’asymptotically 
equal” to the one obtained in the continuous case. The above results are merged in the following 
theorems: 


Theorem 1. Assume A > 0. Then the sequence of minimizers a? 1111 holds 


Ufiim S ^°> b Qj A, strongly in L p (fl) for any p ^ 1, 
u ™m - T - °> b Qj A, uniformly on fl, 


where u/f' n is a locally affine interpolation of «™ ln defined in Lemma 


16 


(14) 

(15) 
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Theorem 2. For A > 0, any sequence of monotone vectors xa satisfying dissipates the 
entropies H a ^ and F QjA at least exponential, i.e. 

H«,a(x1) - H™" ^ (H QjA (x° a ) - H™) e-^^, and (16) 

F q ,a(K) - F™ ^ (F a , A (5&) - F™) e~^ tn , (17) 

mt/iH™ = H QiA (x^ ln ) and F™ A = F a)A (x^ lln ). The associated sequence of densities u a further 
holds 

IK - K in |l!i(fi) < c q , a (H Q , A (xi) - H“) A n ) e^S ( 18 ) 

for any time step n = 1,2,..., where c a ^\ > 0 depends only on a, A. 


Let us now consider the zero-confinement case A = 0. In the continuous setting, the long¬ 
time behaviour of solutions to |l]) with A = 0 can be studied by a rescaling of solutions to ([!]) 
with A > 0. We are able to translate this methode into the discrete case and derive a discrete 
counterpart of [IS Corollary 5.5], which describes the intermediate asymptotics of solutions that 
approach self-similar Barenblatt profiles as t —* oo. 


Theorem 3. Assume A = 0 and take a sequence of monotone 
exists a constand c a > 0 depending only on a, such that 


satisfying (11). 


Then there 


ll w A 


:> A,a,ollL 1 (Q) 


Cr 


H 


«,o(4) 


iJinin 

- n e *,0 


(RaY 


with := (l + a T (2a + 3 )t n ) 


b T (2a + 3) 


where b A a 0 is a rescaled discrete Barenblatt profile and a T ,b T > 0, such that a T ,b T —> 1 for 
t —* 0, see section\3J\for more details. 


Before we come to the analytical part of this paper, we want to point out the following: 
The ideas for the proofs of Theorem [2] and [3] are mainly guided by the techniques developd in 
[45] - The remarkable observation of this work is the fascinating structure preservation of our 
discretization, which allows us to adapt nearly any calculation from the continuous theory for 
the discrete setting. 


1.4. Structure of paper. In the following section [2] we point out some of the main structural 
features of equation ([I]) and the functionals TL a ,\ and F a _\, and show that our scheme rises 
from a discrete L 2 -Wasserstein gradient flow, so that many properties of the continuous flow 
are inherited. Section [3] treats the analysis of discrete equilibria in case of positive confinement 
A > 0: we prove convergence of discrete stationary states to Barenblatt-profiles or Gaussians, 
respectivelly, and analyse the asymptotics of discrete solutions for A = 0. Finally, some numerical 
experiments are presented in section [4] 


2. Structural properties — continuous vs. discrete case 


2.1. Structural properties of equation (jT|) . The family of fourth order equations ([l]) carries a 
bunch of remarkable structural properies. The most fundamental one is the conservation of mass, 
i.e. t b-> ||u(f, OlliAffi) is a constant function for t e [0, +oo) and attains the value M := ||u°||Li(f2)- 
This is a naturally given property, if one interprets solutions to |l]) as gradient flows in the 
potential landscape of the perturbed information functional 


*Fa, A (l^) 



(19) 


equipped with the L 2 -Wasserstein metric W 2 ■ As an immediate consequence, J - a ,a is a Lyapunov 
functional, and one can find infinitely many other (formal) Lyapunov functionals at least for 
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special choices of a — see for a = ^ or mmm for a = 1. Apart from T a , a, one of 

the most important such Lyapunov functionals is given by the A Q ,A-convex entropy 


\{u) = (p a (u)dx + 

Jn 



ip a (s) 


a—1/2 > 

0i/ 2 sln(s), 


ae (|, 1] 
^ _ 1 


( 20 ) 


It turns out that the functionals F a ,\ and H a ,a are not just Lyapunov functionals, but share 
numerous remarkable similiarities. One can indeed see ([Tj) as an higher order extension of the 
second order porous media/heat equation [36] 


8 s v = - gradyy 2 H a ,\(v) = -Q a d xx (v a ) + A a ,\{xu) x , 


( 21 ) 


which is nothing less than the L 2 -Wasserstein gradient flow of H a , Furthermore, the unper¬ 
turbed functional T a , o, be. A = A ai A = 0, equals the dissipation of Hap along its own gradient 
flow, 

•Fa,o(v(s)) = - 4 -n a , 0 {v(a)). (22) 

as 

In view of of the gradient flow structure, this relation makes equation (jTJ) the “big brother” 
of the porous media/heat equation plj ), see [23] 145! for structural consequences. Another as¬ 
tonishing common feature is the correlation of 8F a ,x and H a ,\ by the so-called fundamental 
entropy-information relation: For any u e V(Q) with H a ^\(u) < oo, it holds 

F a ,\(u) = |grad Wa H a ,\\ 2 + (2a - l)A a ,\H a ,\(u), for any A ^ 0, (23) 


see 


, Corollary 2.3]. This equation is a crucial tool for the analysis of equilibria of both 
functionals and the corresponding long-time behaviour of solutions to |l]) and (21). 

In addition to the above listing, a typical property of diffusion processes like ([I) or (21) with 
positive confinement A, A a ^ > 0 is the convergence towards unique stationary solutions zt 00 and 
u 00 , respectively, independent of the choice of initial data. It is maybe one of the most surprising 
facts, that both equations ([T]) and (21) share the same steady state, i.e. the stationary solutions 
it 00 and v°° are identical. Those stationary states are solutions of the elliptic equations 


“( P «( u ))xi + A o,a {xu) x = 0, (24) 

with P Q (s) := 0 a s a+1 ' 2 , and have the form of Barenblatt profils or Gaussians, respectively, see 
definition (12 )&( 13). This was first observed by Denzler and McCann in [23], and further studied 
in |45j using the Wasserstein gradient flow structure of both equations and their remarkable 
relation via (22). 

In case of a e {i, 1}, the mathematical literature is full of numerous results, which is because 


of the physical importance of ([I]) in those limiting cases. 


2.1.1. DLSS euqation. As already mentioned at the very beginning, the DLSS equation - ([I]) 
with a = ^ — rises from the Toom model J2H G5] in one spatial dimension on the half-line 
[0,+oo), and was used to describe interface fluctuations, therein. Moreover, the DLSS equation 
also finds application in semi-conductor physics, namely as a simplified model (low-temperature, 
field-free) for a quantum drift diffusion system for electron densities, see [35] , 

From the analytical point of view, a big variety of results in different settings has been 
developed in the last view decades. For results on existence and uniqueness, we refer f.i. 
0 [28l [35] [301 ESI [39] , and [121 HU HH [301 [381 HU [45] for qualitative and quantitative descriptions 
of the long-time behaviour. The main reason, which makes the research on this topic so non¬ 
trivial, is a lack of comparison/maximum principles as in the theory of second order equations 
plj ). And, unfortunatelly, the abscents of such analytical tools is not neglectable, as the work 
[7] of Bleher et.al shows: as soon as a solution u of ([l]) with a = \ is strictly positive, one can 
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show that it is even C^-smooth, but there are no regularity results available from the moment 
when u touches zero. The problem of strictly positivity of such solutions seems to be a difficult 
task, since it is still open. This is why alternative theories for non-negative weak solutions have 
more and more become matters of great interest, as f.i. an approach based on entropy methodes 
developed in [SUES]. 

2.1.2. Thin-film equation. The thin-ffim equation — ([I]) with a = 1 - is of similar physically 
importance as the DLSS equation, since it gives a dimension-reduced description of the free- 
surface problem with the Navier-Stokes equation in the case of laminar flow, [IS]. In case of 
linear mobility — which is exactly the case in our situation — the thin-film equation can also be 
used to describe the pinching of thin necks in a Hele-Shaw cell in one spatial dimension, and thus 
plays an extraordinary role in physical applications. To this topic, the literature provides some 
interesting results in the framework of entropy methods, see 0 [T8| {29]J. In the (more generel) 
case of non-negative mobility functions m, i.e. 

d t u = — div(m(M) D A u), (25) 

one of the first achievements to this topic available in the mathematical literature was done 
by Bernis and Friedman (4]. The same equation is observed in [5], treating a vast number of 
results to numerous mobility functions of physical meaning. There are several other references 
in this direction, f.i. Griin et. al J5J1221SI3, concerning long-time behaviour of solutions and the 
non-trivial question of spreading behaviour of the support. 

2.2. Structure-preservation of the numerical scheme. In this section, we try to get a 
better intuition of the scheme in section o Foremost we will derive (pi) as a discrete system 
of Euler-Lagrange equations of a variational problem that rises from a Ld-Wasserstein gradient 
flow restricted on a discrete submanifold Vg(SY) of the space of probability measures ’P(Sl) on id. 
This is why the numerical scheme holds several discrete analogues of the results discussed in the 
previous section. As the following section shows, some of the inherited properties are obtained by 
construction (f.i. preservation of mass and dissipation of the entropy), where others are caused 
by the underlying dicsrete gradient flow structure and a smart choice of a discrete L 2 -Wasserstein 
distance. Moreover it is possible to prove that the entropy and the information functional share 
the same minimizer x™ 111 even in the discrete case, and solutions of the discrete gradient flow 
converges with an exponential rate to this stationary state. The prove of this observation is more 
sophisticated, that is why we dedicate an own section (section [3]) to the treatment of this special 
property. 

2.2.1. Ansatz space and discrete entropy/information functionals. The entropies a and Fa,\ 
as defined in (|20|)fc([2|) are non-negative functionals on V(/l). If we first consider the zero- 
confinement case A = 0, one can derive in analogy to [46j the discretization in © of o just by 
restriction to a finite-dimensional submanifold Vg(Cl) of For fixed K e N, the set Vg(fl) 

consists of all local constant density functions u = u^[x] (remember definition ([8])), such that 
x e R^ 1 is a monotone vector, i.e. 

xep := {(io, ■ • •, x K ) \ x 0 < xx < ... < x K -i < x K } £ R K+1 . 

Such density functions u = u^[x] e Vg(ST) bear a one-to-one relation to their Lagrangians or 
Lagrangian maps , which are defefined on the mass grid [0, M\ with uniform decomposition (0 = 
£o , ■ ■. ,£k, ■ ■ ■ ,£k = M). More precicely, we define for the local affine and monotonically 

increasing function X = X^[x] : [0, M ] —» f2, such that X(^) = Xk for any k = 0, ..., K. It then 
holds »oX = ^ for u e 'Pg(il) and its corresponding Lagrangian map. For later analysis, we 
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introduce in addition to the decomposition {£fc}^ 0 the intermediate values (£ fc _i,£s,... ,£ k _i) 
by = §(£fe +£fc-i) for k = 

In view of the entropy’s discretization, this implies using a change of variables 

x = X^[x], and the definition ([7| of the x-dependent vectors z 

r * 

H a ,o(x) = H a ,o(uj[x]) = v? q (u 5 [x]) dx = 6 V f a (z k _ i), 

fc= i 


which is perfectly compatible with <©■ Obviously, one cannot derive the discrete information 
functional F aj o in the same way, since J~„s) is not defined on V&(££). So instead of restriction, 
we mimic property (22) that is for any xe 

F q , 0 (x) = r 1 4H ai o(x) T 4H aj0 (x) = (V5H Q! o(x), V5 H Qj o(x)) 5 ■ 


Here, the fcth component of 5j/(x) holds [<?x/(x)] fc = d Xk f(l) for any k = 0,..., K and arbitrary 
function /:£$—» R. Moreover, we set Va/(x) = (5 _1 <?5f/(x) and introduce for v, w e M A + 1 the 
scalar product (•, -} g by 


K 

<v, w) s = 5 Y, v k w k , with induced norm ||v|| 5 
k =o 



Example 4. Each component z K of z = z^fx] is a function on £ 5 , and 


d*z K = —z: 


2 e K+\ 


(26) 


where we denote for k = 0,..., K by e k e K A+1 the (k + 1 )th canonical unit vector. 


Remark 5. One of the most fundamental properties of the L 2 -Wasserstein metric W 2 on 'P(fl) 
in one space dimension is its excplicit representation in terms of Lagrangian coordinates. We 
refer DEO] for a comprehensive introduction to the topic. This enables to prove the existence of 
K-independent constants Ci,C 2 > 0, such that 


c i l|x- yla W 2 (u 5 [x],U5[y]) ^ c 2 ||x-y|| 5 , for all x, y e pa. (27) 

A proof of this statement for = [a, 6] <z (— 00 , +ao) is given in [36] Lemma 7], and can be easily 
recomposed for = K. 


Let us further introduce the sets of (semi)-indizes 

4 = {0,1,...,*}, and I k 2 = {\\>"->K-\}. 

The calculation (26) in the above example yields the expizit representation of the gradient 


dxH Q , 0 (x) = 0 q 6 ^ 


1/2 


(28) 


and further of the discretized information functional 


z k+i ~ z k-i 


F a ,o(x) = ||V 5 H q ,o(x)||^ = e 2 J Y 

fee 1“. 

In the case of positive confinemend A > 0, we note that the drift potential u i —* |x| 2 u(x) dx 

holds an equivalent representation in terms of Lagrangian coordinates, that is namely X i—> 
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^ |X(£)| 2 d£. In our setting, the simplest discretization of this functional is hence by summing- 
up over all values x k weighted with S. This yields 


H q ,a(x) = H a0 (x) + Y 

fee 1 ° 


and F a , A (x) = F Q>0 (x) + -5 Y 


Ffc 


fcel°„ 


as an extension to the case of positive A, which is nothing else than Note in addition, 

that SJ] keI o K \x k \ 2 = |[x||^. 

A first structural property of the above simple discretization is convecity retention from the 
continuous to the discrete setting: 


Lemma 6. The functional x >-» H QiA is A at \-convex, i.e. 


H a , A ((l - s)x+ sy) < (1 - s)H aiA (x) + sH aiA (y) - ^|A(1 - s)s ||x-y||j (29) 

for any x, y e and s e (0,1). It therefore admits a unique minimizer x^ 1111 e If we further 
assume A q A > 0, then it holds for any ice fj 

A a , A 


||x- xflg < H QjA (x) - H a , A (xf n ) < — I I V * H «,a(x X " 2 


a,A 


<5 ' 


(30) 


Proof. If we prove (29), then the existence of the unique minimizer is a consequence [461 Propo- 
sition 10]. By definition and a change of variables, we get for a e (|, 1] 


H a , 0 (x) = K 


r M f(A sl/2 ne(- ll 

t ,o(u 5 [x]) = V>a(Xa[x] 4 ) d£, with if a (s) = l “ a_1/2 ’ , 

Jo (-Hi/ 2 m(s), a =2 


hence x i —> H Q| o(x) is convex. Since the functional x i—> ||x|b holds trivially 

||(1 - s)x + syf s «S (1 - s) ||x||j + s ||y|| 2 - (1 - s)s||x-y||* 

_^ ^ ^ ^ 2 

for any icjefj and s e (0,1), the functionals H a A (ic) = H aj o(x) + —2 ||xy hold ([29]) . 

Deviding (29) by s > 0 and passing to the limit as s J, 0 yields 

H q , a (x) -H a>A (y) 5 s H Q!A (x)(x-y) - ||x-y||^. 


The second inequality of (30) easily follows from Young’s inequality |o6| < £|a| 2 + (2£) _1 ||6| 2 
with e = (2<fA ajA ) _1 , and even holds for arbitrary ye pj. 

To get the first inequaltiy of (30), we set x = x]f ln and again devide (29) by s > 0, then 

H Q , A ((1- S )^ + S y)-H Q , A (^) A QjA ^ ,2 

--- < LL ajA (y) - tL a , A (x 5 )-— (1 - s) ||x 5 - , 

where the left hand side is obviously non-negative for any s > 0. Since s > 0 was arbitrary, the 
statement is proven. □ 

As a further conclusion of our natural discretization, we get a discrete fundamental entropy- 
information relation analogously to the continuous case (23). 

Corollary 7. For any A ^ 0, every xe y s with H q q(x) < oo we have 

F a , A (x) = ||VH QiA (x)! 2 + ( 2a - 1 ) A a, A H ctjA (x), for a e (i, 1] and (31) 

(32) 


f i/ 2 , A (x) = |V 5 H 1/2iA (x)|| + A 1/2jA , for a=\ 
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Remark 8. Note that the above seemingly appearing discontinuity at a = f is not real. For 
a>\, the second term in right hand side of (311 is explicitly given by 

a- 1/2 . 

(2a - l)A a , A H a , A (x) = (2a - 1 )A q , a I 0 a 5 £ ||x||* 


keT 


1/2 


1/2 ' 2 

= 2A QiA 0 a (5 ^ z k~^ 2 + fta — 1) J —— n^iij , 


. A a A 2 


KSl 


1/2 


For a J, one gets A QjA —» A^^, 0 a —» 5 and especially z k 1 ~ —* M = 1, 77ie 

drift-term vanishes since ( 2 a — 1 ) —> 0 . 

Proof of Corollary Vft ■ Let us first assume a 6 (|,1]. A straight-forward calculation using the 
definition of ||.|| 5 , Ns and djH ajA in (28) yields 

||V s H a , A (*)||= = S- 1 (d ii U a ,x^),ds t ll a ,x(5iy) 

= ||ViH Q , o (x)||^-20 a A Q , A (5 ^ z F~ k + A l,x 6 Ti l**l 2 - (33) 


feel 


1/2 


k€l° 


Here we used the explicit representation of d 5 jH ajA (x), remember 


d 2 H QjA (x) = 0 Q <5 ^ 


z K 


/■cel 


1/2 


fee 1° 


and especially the definition of ([7|, which yields 


W 0 a d ^ 


Zk, 


feel 


1/2 


"i A a? X k G k / 0aA a A (5 Zk 


= -0 a A ajA <5 J] 


a— 7 
' 


/tel 


1/2 


Since a # L^ we can wr ite 20 Q = (2a — 1) a e ^ 2 ■ Further note that the relation A QjA = 
A /A/(2a + 1) yields 


A 2 - 
7 V,a — 


A A / 1 


2a -j- 1 2 \a + 1/2 


1 - 


a — 1/2 A A / 2a- 1 


: + 1/2 


1 - 


2a + 1 


Using this information and the definition of N a ,o> we proceed in the above calculations by 


l|V*H a , A (x)|| 2 = F a , 0 (x) - (2a - l)A a , A H a>0 (x) + \ (l - 6 £ 

' ' k€ 1° 


\Xk\ 


A A 2 

= F Qi0 (x) - (2a - 1 )A QjA H Qi0 (x) + -5 |a; fc | 2 - (2a - 1)—^<5 ^ \x k \ 


fee 1 ° 


fee 1° 


= F a , A (x) - (2a - l)A aiA H ajA (x). 


In case of a = we see that © 1/2 = and A]y 2)A = 
|V 5 H 1/2jA (x)g = |V 5 H 1/2>0 (x)|g - A 1/2jA <5 ^ 

K ei^ 2 


y/A/2. We hence conclude in (33) 

z l + S l^l 2 = F «a(x) - Ai /2)A 

fcel° 


□ 
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For the following reason, the above representation of F QjA is indeed a little miracle: From 
a naive point of view, one would ideally hope to gain a discrete counterpart of the funda¬ 
mental entropy-information relation (23), if one takes the one-to-one discretization of the L 2 - 
Wasserstein metric, which is (in the language of Lagrangian vectors) realized by the norm 
x >—> W 2 (u^[x], [x]) instead of our simpler choice x >-* ||x|| 5 . Indeed, with this ansatz, the above 
proof would fail in the moment in which one tries to calculate the scalar product of e^H^o and 
&VV 2 (uj[x], u 5 [x]). This is why our discretization of the L 2 -Wasserstein metric by the norm ||- || tf 
seems to be the right choice, if one is interested in a structure-preserving discretization. 

Corollary 9. The unique minimizier x ^ 1111 e ya of H QjA is a minimizer of F ajA and it holds for 
any x e y s 


Fq.,a(x) - F Q , A (xf n ) < ^±1 ||V S H Q , A (*)||* . 


(34) 


Proof. Equality (31) and 2 q — 1 ^ 0 shows, that x h-» F QjA (x) is minimal, iff ||V 5 H ctjA (x )|| (5 = 0 
and H q A (x) is minimal, which is the case for x = x^ 1111 . The representaion in (31) further implies 

F q>a (x) - F a , A (xf n ) 

= ||V*H aiA (x)||2 - ||V 5 H Q , A (xf 1 )|g + (2a - l)A Q>A (H a , A (x) - H Q , A (xf n )) 


= |V 4 H a , A (x)||* + (2a - l)A a , A (H Q>A (x) - H QjA (xf n )) < 
where we used (|30|) in the last step. 


1 + 


2a- 1 


||V d -H ajA (x) 


I <5 > 


□ 


2.2.2. Inerpretation of the scheme as discrete Wasserstein gradient flow. Starting from the dis¬ 
cretized perturbed information functional F QiA we approximate the spatially discrete gradient 
flow equation 


1 = —V5F QjA (x) 


(35) 


also in time, using minimizing movements. To this end, remember the temporal decomoposition 
of [0, + 00 ) by 

{0 = t 0 < h < ■ ■ ■ < t n < ■ ■ •} , where t n = i„_ 1 + r„, 

using time step sizes t := {tl,T 2, ... ,r„,...} with r n < r and r > 0. As before in the in¬ 
troduction, we combine the spatial and temporal mesh widths in a single discretization pa¬ 
rameter A = (t;i5). For each yeft, introduce the Yosida-regularized information functional 
F a (-,-,-,y) : [0, +oo) x [0, t\ x y 5 by 


F a (A, <7, x, y) = ||x — yII 5 + F ajA (x). 


(36) 


A fully discrete approximation (x ^)” =0 of (35) is now defined inductively from a given initial 
datum x^ by choosing each x^ as a global minimizer of F Q (A, r„, •, x^ -1 ). Below, we prove that 
such a minimizer always exists (see Lemma 11). 

In practice, one wishes to define x^ as - - preferably unique — solution of the Euler-Lagrange 
equations associated to F a (A, r n , •, x^ 1 ), which leads to the implicit Euler time stepping: 


V n-l 
X A 


= —V,yF a)A (x). 


(37) 


Using the explicit representation of it is immediately seen that (37) is indeed the same 

as <©■ Equivalence of ( |37[ ) and the minimization problem is guaranteed at least for sufficiently 
small r > 0, as the following Proposition shows. 
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Proposition 10. For each discretization A and every initial condition x° G £s, the sequence of 
equations (371 can be solved inductively. Moreover, if r > 0 is sufficiently small with respect to 6 
and F ay \(jp), then each equation (37) possesses a unique solution with F^^x) < F QiA (x°), and 
that solution is the unique global minimizer of F Q (A, r n , ■, x^~ ). 


The proof of this proposition is a consequence of the following rather technical lemma. 


Lemma 11. Fix a spatial discretization parameter 6 and a bound C > 0. Then for every y G £s 
with F a , A (y) ^ C, the following are true: 

• for each a > 0, the function F Q (A, a, •, y) possesses at least one global minimizer x* e y^; 

• there exists arc > 0 independent ofy such that for each a e ( 0 , tc), the global minimizer 
x* G is strict and unique, and it is the only critical point of F Q (A, a, •, y) with F ajA (x) < 

C. 


Proof. Fix y e fj with Fo, jA (y) < C , and define the nonempty (since it contains y) sublevel 
Ac := (F a (A, er, -,y)) ([0, C]) <z y^. First observe, that any xg Ac holds 

V2 aC 3* ||y — x || 5 3* <^||y-x||oo 3* Hyloo - ||x||oo|, (38) 

hence Ijxljoo is bounded from above by f2aC + ||y||oo- Especially, 

max|xfc| < V2<5- 1 ctC + ||y||ao =: L(8,cr,f), (39) 

ke I° K 

which means, in the sense of density functions, that any u = 115 [x] with x G Ac is compactly 
supported in [— L{5, a, y), L(S, a, y)]. Consequently, take again x G yj arbitrarily and declare 
z* = min KeI i /2 z K and z* = max KgI i /2 z K , then on the one hand the conservation of mass yields 
the boundedness of z * from above, 

1= f u 5 [x]dx= 2 z~ 1 (x K+ i-x K _i)^2L(8,a,y)(z*)- 1 , 


and on the other hand F ajA (x) < F QjA (y) < C yields an upper bound for z*, as the following 
calculation shows: 


(**r»-(*.r»« E ( S*| (*S 

hB'£ 

/ . 1 \ l/(a + l/2) 

** < (M0- 1 C+(2L(<5,a,y)) Q+5 ) 


Z t — Z f 

/€+ 4 K— 


(40) 


Collecting the above observations, we first conclude that Ac c y s is a compact subset of R A+1 , 
due to |xo|, \xk\ ^ L(S, a, y) and the continuity of F Q|A . Moreover, every vector x G Ac satisfies 
x K+ i — x K _i (g <5(z*) _1 Jg x for all k g with a positive constant x that depends on C and 
L(5,a,y). Thus Ac does not touch the boundary (in the ambient K A+1 ) of y^. Consequently, 
Ac is closed and bounded in y^, endowed with the trace topology. 

The restriction of the continuous function F a (A, a, -,y) to the compact set Ac possesses a 
minimizer x* G Ac- We clearly have F QjA (x*) < F a>A (y) C, and so x* lies in the interior of 
Ac and therefore is a global minimizer of F a (A, a, -,y). This proves the first claim. 

Since F Q A : y$ —» M is smooth, its restriction to Ac is Ac-convex with some Ac < 0, i.e., 
A (x) > Aclif+i for all x G Ac- Independently of y, we have that 

£ 

5|F a (A,o-,x,y) = 5|F a A (x) + - t K+1 , 

T 
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which means that x <—> F Q (A, cr, x, y) is strictly convex on Ac if 

S 

0 < a < tq := 


(-Ac) 

Consequently, each such F Q (A, cr, •, y) has at most one critical point x* in the interior of Ac , and 
this x* is necessarily a strict global minimizer. □ 


Remark 12 (propagation of the support). Take a solution x A of (36) with density functions 
ua- As we already noted in the above proof, any density u\ has compact support in [— L n ,L n ] 
with L n = L(5,r n ,x^ -1 ) as in (39). Hence 


'■A II00 


^25- + K-'Hoo 

||xa||oo < V 2 S - 1 ■\J T jF a ,x{^A 1 ) + I^aIIoo, 


J = 1 


which is the best we can assume in case of A = 0. If A > 0, one can find a much better bound on 
the support of ua, namely by replacing (38) by, 

2S- 1 


IIxaIIoo ^ ||x£||j 


< 


2A- 1 

F«,a(xa _ 1 ) < -^F a , A (x° A ). 


A A 

2.3. Some discrete variational theory. In this section, we consider an arbitrary function 
V : y^ —* (—00,+00] and assume the existence of a value r* = r*(V) > 0, such that for any 
xe the minimization problem 


v er(y, x) := ||y — x|j + V(y) 


(41) 


has a solution x a in for any ex e (0, t*]. In the literature, the function cr x CT is known as 
the De Giorgi’s variational interpolant connecting x and x T *, see f.i. dj section 3.1]. Another 
interesting object in this context is the discrete Mureau-Yosida approximation a <-> V cr (x cr ,x), 
which can be defined for any xe f^. Analogously to the theory developed in pQ section 3.1], one 
can even introduce a discrete version of the local slope of V, i.e. 


| <5,5 V | (x) := limsup 

y€jj:y->x 


(V(x)-V(y))' 


(42) 


l|x-y|li 

It turns out that V| (x) = || VsV(x) || 5 , which is a consequence of the lemma below and an 
analogoue calculation done in PQ Lemma 3.1.5]. 

The above definitions remind of their continuous counterparts as defined in [I, section 3.1], 
and the reader familiar with pQ knows, that those objects are well-studied. Some properties of 
the discrete Mureau-Yosida approximation, which will be needed in later sections to study the 
asymptotic behaviour of solutions to are listened in the following lemma. The proof is a 
special case of pQ Theorem 3.1.4 and Lemma 3.1.5]. 

Lemma 13. Fixxe y^ and declare by x CT the De Giorgi’s variational interpolant. Then it holds 
for any a e (0, r*] 


2er 


f 


2r 2 


■ dr = V (x) — V (x CT ). 


(43) 


If we further assume the continuity of a > djV(x cr ) and a <-» <5|V(x CT ), and the validity of the 
system Euler Lagrange equations 

1 


:( J 


?) = -V 5 V(x a ), 
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for a e (0, r*] ; then 


||V 5 V(x)|| 2 = lim 


V.IIS 


crjO 


(7 Z <r|0 (7 crjO 

3. Analysis of equilibrium 


V(x) - V(x cr ) 




(44) 


In that which follows, we will analyze the long-time behaviour in the discrete setting and will 
especially prove Theorem [2] As we have already seen in [45] . the scheme’s underlying variational 
structure is essential to get optimal decay rates. Due to our structure-preserving discretization, 
it is even possible to derive analogue, asymptotically equal decay rates for solutions to (36). 

3.1. Entropy dissipation — the case of positive confinement A > 0. In this section, we 
pursue the discrete rate of decay towards discrete equilibria and try to verify the statements in 
Theorem [2] to that effect. That is why we assume henceforth A > 0. 

Lemma 14. A solution xa to the discrete minizing movement scheme © dissipates the en¬ 
tropies H a j and F QiA at least exponential, i.e. 

(1 + 2r„A) (H q , a (x£) - H™ n ) < H Q , A (5T X ) - , and 

(1 + 2t„A) (F a , A (x£) - F™) < Fa,^- 1 ) ^ F 

for any time step n = 1 , 2 ,.... 


pmin 

, A 


(45) 

(46) 


Proof. Due to (31), the gradient of the information functional F a>A is given by 

d 2 F ajA (x) = 2r 1 (d 5 H QjA (x)) T d|H Q , A (2) + (2a- l)A a , A d*H a , A (x), 
which yields in combination with the A ajA -convexity of H Q]A and ( |37[ ) 

H q , a (^ _1 ) - H a , A (x^) 

>T„<V 4 F aiA (x£),V 4 H a , A (x£)> 4 

^2r„ <V 4 H a , A (x£), 5lH a , A (x)V 5 H a , A (x^)) 5 + r„(2a - l)A a , A ||V 4 H a , A (x£)|| 2 

^2r„A ajA ||V 4 H a , A (x£)||* + r n (2a - l)A a , A ||V 4 H Q , A (x^)||^ ^ r„(2a + 1 )A q , a ||V 4 H q , a (x^)||^ . 

(47) 


Using inequality (30), we conclude in 

(1 + 2r„(2a + 1)A* >a ) (H a , A (x£) - H a , A (2f n )) < H*,^" 1 ) - H Q , A (*f n ) 

for any n = l, 2,_ Since (2a+ 1 )A^ >a = A, this shows ( [45] ) . To prove ( [46] ) , we proceed as above. 

First introduce for er > 0 the vector x a as unique minimizier of y >-» A. ||y _ x™ || y H ttiA (y). 
Then the fundamental property of the minimizing movement scheme reads as 


F Q , A (xl) - F«, a (x ct ) < (|| 


-fn —11 


ll-K — x rl_ 1 H 2 ') 

|| X A X A || 4 J 


< 7T- ll x <T - Xa|| 4 (||x ct - 1 \\ s + ||Xa 


-fn — 1II 


Devide both side by a and pass to the limit as a J, 0, we get 

<V 4 F a>A (x£), V 4 H a , A (x£)> s < - ||V 4 H q , a (xa)| (5 ||xa - *a“ 

7"t7. 


\\S ’ 


due to (44), and thanks to (47) further 


-*n — 11| 


r n (2a + 1)A QjA ||V 4 H a , A (x£)|£ < ||V 4 H aiA (x£)||, ||^ - 


II <5 ' 
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-?" -1 II 

X A IU 
r+n— 1 ||2 


and 


(48) 


As a consequence, we get two types of inequalities, namely 

W(2a+l)A ||V«H a , A (S£)|| 4 < ||- 
2 t 2 A(F q>a (* 2J - F a , A (xf n )) < ||S£ - || 5 , 

where we used A a>A = ^/A/(2a + 1) and (f34|) . To get the desired estimate, remember the De 
Giorgi’s variational interpolation: Fix x^ -1 and denote by x^ a minimizer of y >—* ^ ||y — |“+ 

F a , A (y) for a e (0 ,t„]. Then xJJ connects x^ _1 and and the monotonicity of a F ajA (x^) 
and (48) yields for any a e (0, r] 

2cr 2 A(F Q A (x^) — F a A (x“ in )) ^ 2r 2 A(F a A (x^) — F a A (x“ in )) 

i,2 (49) 


< I Is?” _ •? r!_1 l, 

| X er X A ||<5 • 


Further remember the validity of (|43|) in Lemma 13 which gives us in this special case 

f 

Jo 


Fa,A(*£) + 


II-? 71 — y 71-1 !! 

-A-A -A- A 


Hx^-xT 1 ! 


2 T n 


+ 


2a 2 


■da = F aiA (xST 1 ). 


Inserting (J49J) in the above equation then finally yields 

(1 + 2r n A) (F q , a ( 5%) - F a , A (xf n )) ^ F QjA (x^ _ 1 ) - F«, A (xf n ), 
and the claim is proven. 


□ 


Remark 15. In the continuous situation, the analogue proofs of (45) and (46) require a more 


deeper understanding of variational techniques. An essential tool in this context is the flow inter¬ 
change lemma, see f.i. |45( Theorem 3.2]. Although one can easily proof a discrete counterpart 
of the flow interchange lemma, it is not essential in the above proof, since the smoothness of 
x i—» H q A (x) allow an explict calculation of its gradient and hessian. 

Lemma [l4] paves the way for the exponential decay rates of Theorem [2j Evectively, (16) and 
( |17| ) are just applications of the following version of the dicrete Gronwall lemma: Assume {c n } ne N 
and {y n }n €n to be sequences with values in R + , satisfying (1 + c n )y n y n -i for any n e N, then 


_y»-i c k 

y n yoe ^ k=0 1+Ck , 


for any n e N. 


This statement can be easily proven by induction. Furthermore, inequality (18) is then a corollary 


of (16) and a Csiszar-Kullback inequality, see tn, Theorem 30]. 


3.1.1. Convergence towards Barenblatt profiles and Gaussians. Assume again A > 0. It is an 
striking fact (see [23]), that the stationary solutions u 00 and u 00 of |l]) and (21), respecively, are 
identical. Those stationary states have the form of Barenblatt profils or Gaussians, respectively, 


b Q ,A = (a - b\x\ 2 ) 
b 


2 \ 1/(a —1/2) , a — 1/2 

' , b = — — A a A if a > 1/2 and 

s/2a 


= ae Ai / 2 aI x I if a = 1/2, 


A/ 2 ,a — o.e 

where a e R is chosen to conserve unit mass. 

To prove the statement of Theorem [l] we are going to show that the sequence of functionals 
fL s a A (-) : V(n) — (-oo, +oc] given by 


<a(«) == 


\Ha,\(u) ,ueVg(Cl) 


+ 00 


, else 


T-congerves towards Ha, a- More detailed, for any u e 'P(fl) it holds 

(i) liminf5_ > o ^i, A ( Ui ) ^ ^a,A (u) for any sequence us with lim,5_>o W 2 (u$, u) = 0. 
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(ii) There exists a recovery sequence us of u, i.e. limsup^Q TL s a x (us) < TL a ,\(u) and 
lim^_>o W 2 (u $, u) = 0 . 

The T-convergence of TL 5 a A (-) towards TL a ,x is a powerful property, since it implies convergence 
of the sequence of minimizers w“ ln = u^fx^ 1111 ] towards b QiA or b]/ 2 iA , repsectively, w.r.t. the 
L 2 -Wasserstein metric, see [ 8 ]. To conclude even strong convergence of u™ ln at least in L P (Q) 
for arbitrary p ^ 1, we proceed similar as in [461 Proposition 18]. Necessary for that, recall that 
the total variation of a function / G T 1 (S1) is given by 


{/}tv := sup | f(x)ip'(x) dx 


ip e Lip(fl) with compact support, sup || < 1 

xefl 


(50) 


we refer m Definition 1.1]. If / is a piecewise constant function with compact support [x’o, Xk\ , 
taking values f k _ i on intervals ( x k -i,Xk ], then the integral in (50) amounts to 


~ K K -1 

f(x)ip\x) da; = 2 [f(x)ip( ®)]“l ~°_ 1+0 = J] (fk -§ - /k+iM^fc) + o) - f K _iip(x K )- 

k=1 k= 1 


Consequently, for such /, the supremum in (50) equals 

K -1 


{/}tv = Xj l/fc+i - A-i| + l/|l + \fx-i\ 


(51) 


k= 1 


and is attained for every ip e Lip(D) with p(x k ) = sgn(f k — f k +i ) at k = 1,..., — 1, p(xf) = 

sgn(/i) and ip{x K ) = -sgn(f K _i). 


Lemma 16. For any a G [ '(, 1], assume x ^ 1111 e to be the unique minimizer o/H Q A and declare 
the sequence of functions [x]f ln ]. Then 

uf 11 b QjA , strongly in L p (n) for any p > 1, (52) 

—°> b QjA , uniformly on f2, (53) 


w/iere w™ 


n 


1/2 


a local affine interpolation of on f 1, such that for any k elf. 


u4 


(uf n oX 5 [xf n ])fe) = ^ 


\z K , 

\ 2(^+5 "b 5)) 


KG if 

«el§f 


Proof. We will first prove the T-convergence of TL'f A (-) towards TL a ,x- The first requirement (i) 
is a trivial conclusion of the the lower semi-continuity of TL a ,x- For the second point (ii), we fix 
u G V(fl) and assume X : [0, M ] — > [—00, +00] to be the Lagrangian map of u. Further introduce 
the projection map ns : X —» £5 on the space of Lagrangian maps, 

K 

tt 5 [X] = ^ X(&)6>*(0, 

k =0 


where d K : [0, M] —» D are local affine hat-functions with 9 k (lS) = S L K for any k,(G if u I^. We 
claim that the corresponding sequence us to 7 r^[X] is the right choice for the recovery sequence. 
To prove the convergence in the L 2 -Wasserstein metric, we fix e > 0 and take a compact set 
K. ci \—L, L\ c Q with f |x| 2 rt(x) dx < e and \ /c \x\ 2 us(x) dx < e, which can be done due to 
the boundedness of the second momentum. Since X and 7 r^[X] are monotonically increasing, it 
holds for any £ G [0, M\ 

|X(£) - 7rs[xm\ < (X(£ fc ) - X(a-i)), with £ g [£ fc _i, £ fc ], k = 1,..., K 
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and further for S ^ e(2 L) 2 

||X - ^[X]Il 2 ([ 0 ,m]) 

<|X — 7r <5[X]||i2(x-i([-L,L])) + ||X - 77,5 [X] Il2([ 0 jv,/]\x- 1 ([-L i L])) 

^2L5 I |X(4)-X(4-i)| 

X« fc )e[-L,L] 

+ (ll X IU 2 ([0,M]\X-i([-L,L])) + l 7r <5[X]||i2([ 0l M]\X-i([-L,I,]))) 
<4 L 2 S + 8e^ 9e. 


This shows us —> u in L 2 -Wasserstein, as 5 —» 0. The second point (2) easily follows by using 
Jensen’s inequality, 



Taking the limes superior on both sides proves limsup 5 ^ 0 \( u s) ^ 'Ha,xW) an d since 'H ol ^\ 
is lower semi-continuous, we especially obtain lima^o'H^ A (ua) = H ajA (u). Due to the the equi- 
coercivity of 'H s a A (-) and vad U EV{o.) 'H s a x (u) = H ajA (x^ lm ), w™ converges towards b aj > in w.r.t. 
W 2 , by [3 Theorem 1.21] 

The convergence of H Q!A (x™ ln ) to T~L a ^\{ b QiA ) yields on the one hand the uniform boundedness 
of H Q , jA (i^ lln ) w.r.t. the spatial discretization parameter <5, and on the other hand the uniform 
boundedness of F QjA (x]f ln ), which is a conclusion of © and V, 5 H Qi a(x) = 0. Similar to [46| 
Proposition 18], one can easily prove that the term F Q , A (x^ lln ) is an upper bound on the total 
variation of P a (it™ m ) with Pq,(s) := 0 a s a+1 / 2 : Take any arbitrary ye R^ 1 with ||y || oo < 1 and 
define Y = X^[y]. Then 


<V 5 H a , 0 (5f n ),y> 5 = <V 5 H QiA (xf n ),y> 5 - A q , a <xf n ,y> 5 , (54) 

and the left hand side can be reformulated, using ( |28| ) and a change of variables, i.e. 

Vk+ i Vk —i 


<V 5 H Q , 0 (xf n ),y> 6 = -«5 ^ P a (z K )' 


I 


«el' 

X K 


1/2 


r>M 

1 1 


p Q ( U r n °xs m )d 4 Yde 


P a (uf m )p x dx 


with the Lipschitz-continuous function ip : \xq,Xk\ —» [—1,1] defined by ip{x) := X^[y] o 
(X™ n ) _1 . Since w™ ln is equal to zero in fi\[xo, %k], we can dehne any extension : ST —> [—1,1] 
of (p with compact support, and exchange it in the above calculation without changing the value 
of the integral. Moreover, by the Cauchy-Schwarz inequality, ||y|| 5 < ||y||oo ^ 1 and (541, 


f P a (uf n )^dx ^ |V 5 H a , A (xf n )L + A q , a \\lf 

Jn 
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which is uniformly bounded from above, due to (31) and the uniform boundedness of 

and F Qj ^(xg lm ). This proves the uniform boundedness of {P a (it™ ln )}TV and using the superlinear 

growth of s « P a (s) and m Proposition 1.19], we conclude in |52| . 

To proof (53), we show that the H 1 ( B)-norm of zi™ ln is bouned by the information functional. 
This was already done in m for a = i, where we also showed 

Zu , 1 + Z ; ' / N 2 

||- mm || 2 _ X V '"+5 

\\ u s llffi(n) — 0 2-i 

fce l° K 

So assume a G (g,l], then the concavity of the mapping s i-» s a ~ 1 ^ 2 yields for any values 


k— ^ 


z k+ i 


k— ^ 


b > a > 0 


- (a+\ 


rb 

i) 

J a 


ds ^ (a + i)- 


"2 + a 


-(b-a), 


and further (b a+ 2 — a a+ 2 ) 2 ^ (a + \) 2 - -^-(6 — a) 2 . Therefore 

|| „-min || 2 _x V Zk +h +Zk ~h f Z k+ k ~ Z k- 5 

II U s llHi(n) - 6 2j - 2 

fcel° 


< |U'/ rni n|l 2 ( 1—a ) A V 

ll u <5 lli®(n) 0 2j 
fcel° 


v <5 

2a-l . 2a—1 

“fc+ i + 2 fc- i 






< 


2 || 


|| 2(1—a) 


(a + -) 2 S ^ 
+ 2 1 fc e l° 


Z k+ A 2 fe-i 




4a||: 


|| 2(1—a) 


(a+i) 2 


Fa,A(5f n )- 


□ 


3.2. The case of zero confinement A = 0. We will now consider equation 0 in case of 
vanishing confinement A = 0, hence 

d t u = — [u(u a ^ 1 uf x ) x ) x , for (t, x) G (0, +oo) x B, (55) 

and it(0) = u° for arbitrary initial density u° e From the continuous theory, it is known 

that solutions to (551 or (21) with A a ,\ = 0 branches out over the whole set of real numbers 
S2 = R, hence converges towards zero at a.e. point. This matter of fact makes rigorous analysis of 


the long-time behaviour of solutions to (55) more difficult as in the case of positive confinement. 


However, the unperturbed functionals Haft and F a ,o hold the scaling property, see again [do] . 

W a ,0 0r«) = 7’~ (2a_1)/2 'Ha,o(w), and Ja,o(frw) = r~ {2a+1) T afi (u ), (56) 

for any 7' > 0 and d r u(x) := r~ 1 u(r~ 1 -) with u e 'P(fl). Due to this, it is possible to find weak 
solutions to a rescaled version of (55) by solving problem ([!]) with A = 1. More detailed, it holds 
the following Lemma: 

Lemma 17. A function u G L 2 OC ((0, T); W 2 ’ 2 (tt)) is a weak solution of 0 with A = 1, iff 


w(t, •) = D RW u(log(R(t),-), with R(t) := (l + (2a + 3)f) 


l/(2a+3) 


(57) 


is a weak solution to (55). 


A consequence of the above Lemma is, that one can describe how solutions in to (55) vanishes 
asymptotically as t —» oo, although the gained information is not very strong and useful: In 
fact the first observation (without studying local asymptotics in more detail) is, that w decays 
to zero with the same rate as the rescaled (time-dependent) Barenblatt-proffle b* 0 defined by 
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b* o(f, •) := with R(t) of (57). It therefore exists a constant C > 0 just depending on 

R a ,o{w°) = Ha,o{u°) with 

\Ht,-)-Ko(^-)\\L^n)<CR(t)-\ (58) 

for any t > 0. In [45j . this behaviour was described using weak solutions constructed by min¬ 
imizing movements. We will adopt this methodes to derive a discrete analogue of (581 for our 
discrete solutions xa of (36). 

First of all, we reformulate the scaling operator D r for fixed r > 0 in the setting of monoton- 
ically increasing vectors x e Since d r u(x ) := r _1 u(r _1 -) for arbitrary density in the 

same can be done for ug = u^[x], hence 


K 


„-i; 


K 


T>rU S (x) = V -!(* 1 , x . 1 ](r 1 x) = V 


-I, 


(rx _i,rx , i] 

K O K ' o 


{x) = u 5 [rx](a:) 


^_^x k -x k -i K s - c + 2 j j^rxk-rxk-i 

for any x e ft. The natural extension of D r to the set y^ is hence 

d r x := rx, with corrseponding D r z = z^flVx] = r _1 z. 

As a consequence of this definition, we note that it holds the discrete scaling property for H Qj0 
and F ai o, he. for any r > 0 and xe yj, 


H aj0 (l>rx) = r (2a 1)/2 H Qj0 (x), and F Qi0 (t) r x) = r ( 2 “ +1) Ja,o(x). 


(59) 


The first equality holds, due to H aj0 (x) = H a ,o(uj[x]) and the scaling property (56) of the 
continuous entropy functions. The analogue claim for F Qj0 in (59) follows by inserting d r x into 
d 5 fH Qj o and using d r z = r~ x z, then 


C^xHa,o(^rX) (^r^j 


^±5 = r -(a+V2 )^ Hqi0 (x) 


kgI' 


1/2 


=> F Q ,o(fri?) = ||V 5 H a ,o(3r2)|i5 = r-( 2Q+1 ^ ||V A -H Q , 0 (x) f s = r ~^ 2a+1) F a , 0 (x). 

This scaling properties can now be used to build a bridge between solutions of discrete mini¬ 
mizing movement schemes with A = 0 to those with positive confinement. The following Lemma 
is based on the proof of Theorem |451 Theorem 5.6], but nevertheless, it is an impressive example 
for the powerful structure-preservation of our discrete scheme. 

Lemma 18. Assume x* e y s with F aj >,(x*) < oo. Further fix r > 0 and R> S > 0. Then x e y^ 
is a minimizer of 


y l-> F Q (A,t, y,x*) = ||y ■ 

if and only if dux e $g minimizes the functional 


■S*ll« + Fa,o(y) + ^ IlyllL 


F a (A,f,w,D s x*) = ||w- 0 s x*||^ + F a>0 (w) + ^ llw" 2 

Zr ’ 


with, 


= rSR 2a+2 , A = 


2 ii ■■ ik5 > 

5(1 + At) - I? 


(60) 


(61) 


tR 


It is not difficult so see, that this Lemma holds for all functionals V : fg 


with the same 


scaling property as F a ^ 0 hr (591. 

Proof of Lemma \ 1 A] . To simplify the proof, we show first that we can assume 5=1 without 
loss of generality, which is because of the following calculation: If for R > S > 0 the vector dux 
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minimizes (61), then the linearity of || • || 5 and (59) yields 

S 2 


-lr. *11 2 , c-2a + l-p /o-l* ,-T\ i c2^IIc-1t ->" 2 


F a (\,?,V R x,<) s x*) = — |5- 1 0 ii x-x*|; + S-^F^S-^rx) + S 2 - ||S'- 1 0 Je 3t ||ff 


= g—(2a + l) 


2tS-^+ 3 ) IM- **ls + + S 2a+3 ^ |M|" 


= 5-( 2a+ i) Fa (A, t, 0^x,x*), 
with R = ^ > 1 > 0 and the new constants 

t = TSR 2a+2 S~V a+ V = rR 2a+3 , and A = S 2a+3 + X j} ~ R/S = (1 + At) - £ 

t R/ S tR. 

hence D^x minimizes F a (A, r, d ^x, x*). 

So assume S' = 1 and R > 1 in (61) by now. Further introduce the functional g : x K. —> K 
9(y,r) ■■= ^IKy-x*||j + ^a,o(y) + ^(1 + Ar-r) ly^, 

then by definition 

T~ l g(y, 1) = F a (A,r,y,x*), and (TR 2a+2 )~ 1 g(y,R) = F a (A,r,i) H y,x*). (62) 


For fixed y e y^, a straight-forward calculation shows, that the derivative of r 1 —> g(y,r) holds 


drg(y,r) = (Q r y - X* ,y) s + T a ,o(y) ~ £ ||y||a + ^(1 + Ar^r) |jy||^ 


= -<x*,y> 5 + ^a,o(y) + ^(1 + At) ||y||^ ^ ||y • 

= g{ y- 1 ) ^ ^ II s * II* ■ 


-** || 2 
X* L 


7, + Ja,o(y) + Y ||y||* 


Hence, if x minimizes (601, then the same vector minimizes y <-* g{ y, 1) and further y h-» d r g{ y, r) 
for any 7 ' > 0. By integration 

g(y,r) -g( y,l) = J d s g(y,s)ds = (r-l)(g{y,l) - * ||x*|||) 

=► g(y, r ) = r g(: y» i) - ( r -!) ^ II s * I1* 

for any r > 1 and y e y*. This means especially, that for arbitrary r > 1, the function g(y, r) is 
minimal if and only if g{ y, 1) is so. This proves in combination with (62), that d R x is a minimizer 
of (61). By integration of d s g( y, s) over [r _1 ,1], r > 1, one can analogously prove, that if x e y* 
is a minimizer of (61), the rescaled vector d R -ix has to be a minimizer of (60). □ 

Before we prove the claim of Theorem [3j let us introduce the rescaled discrete Barenblatt- 
profile. Define inductively for n = 0,1,... 


Sj. := 1, £? = (! + r n )S? 


n —1 


(63) 


Further take the minimizer x^ 1111 e y* of the functional x i—> H ail (x). Then denote the scaled 
vector b\ a 0 := and define its corresponding density function b^ a 0 = u*[b^ 0 ], This 

function can be interpreted as a self-similar solution of ( {6T| ) with initial density u^fic^ 1111 ], A = 0 
and with time steps ? n inductively defined by ? n := T n S^ r3 S 2a+2 . 
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Proof of Theorem [5] . As already mentioned above, we define a sequence of functions S'™ induc¬ 
tively trough ( |63| and declare a new partition of the time scale [0, + 00 ) by 

n 

{0 = s 0 < 5 2 < ... < s n < ...}, where s n := ^ r k and f k := r fc S* _1 (S*) 2a+2 , (64) 

k= 1 

and we write t? = (ti, T2,...). As a first conequence of the iterative character of the above 
object, we note that (1 + x) < e x causes S” < e‘" for any n = 0,1,_ Moreover 


Sn = £ T k S k T -\S k ;) 2a+2 = £ r fe (l + r fc ) 2a+2 (^- 1 ) 2 “+ 3 < (1 + r) 2a+2 £ T k e^ 2a+3 ^-P 


k= 1 


fc=l 


fc=l 


This is nice, insofar as the right hand side is a lower sum of the integral (l + r) 2a+2 e (2a+3)s 

hence 


s n (1 + r) 2a+2 (2a + 3) _1 [e^ 2a+3 ) tr * - l] 


(65) 


=► e~ tn < (1 + a T s n {2a + 3)) 1/(2 “ +3) , 

with a T = (1 + T )-( 2a + 2 ) converging to 1 as r —> 0. For a given solution xa of (60) with A = 1 
and fixed discretization A = (r;5), it is a trivial task to check, that the recursively defined 
sequence of vectors is a solution to (61) for S = S™ -1 , R = S™, A = 0 and t = f n defined 

in (64). Henceforth, we write x]~ = Os^^a with the discretization A = (r?; S) We can hence use 
the discrete scaling property of H a j and invoke ( |45| ) of Lemma [14} then 

(1 + 2r„)(S™)^(H a>1 (x|) - H Qil (bX, a , 0 )) < (5 T ™- 1 )% & |H Q , 1 (x| -1 ) - H Qjl (b£+)) 
=> (1 + 2r„)(l + T n )^ (H Qi1 (x|) - H Qil (bl iai0 )) < H Q!l (x^ -1 ) - H Qil (b^ 0 ) 

=> (1 + 2t„)(H Qj1 (x|) - H ail (bi iQi0 )) ^ H^fxl- 1 ) - H ail (b™^ 0 ), 

(66) 

where we used in the last step (1 + r n ) > 1. As before in the proof of (16) of Theorem [2} this 
yields for any n = 0,1,..., due to 


H Qil (x|) - H Qil (b^, Qj0 ) ^ (H Qi1 (x|) - H ail (xf n ))e-a^ 

^ (H ail (x|) - H ajl (xf n )) (1 + a T s n (2a + 3 ))“VW, 


with b r = 1 + 2 r. Theorem |3| follows, using b^ ?a?0 = u«5[bA, a ,o] and a Csiszar-Kullback inequality, 
see da Theorem 30]. □ 


4. Numerical experiments 

4.1. Non-uniform meshes. An equidistant mass grid — as used in the analysis above — leads 
to a good spatial resolution of regions where the value of u° is large, but provides a very poor 
resolution in regions where u° is small. Since we are interested in regions of low density, and 
especially in the evolution of supports, it is natural to use a non-equidistant mass grid with an 
adapted spatial resolution, like the one defined as follows: The mass discretization of [0, M] is 
determined by a vector S= (Co, ■ ■ •,Cat- 1, Cif), with 0 = Co < Ci < • ■ ■ < Cr-i < Cif = M 
and we introduce accordingly the distances (note the convention C-i = Ca+i = 0) 

= C K +i - Ck- 1) and S k = -(£ fc+ j + 5 k _i) = — (C*+i - Cfe-i) 
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for k e l](- 2 and k e 1^-, respectively. The pi 

^ cr u^A-i j g now gj ven ]-,y 

u(x) = z K for x K _i < x < x 


___ .. _ _ .. _ _ A ,_j._ v . _piecewise constant density function u e Vt(i 2) 

corresponding to a vector xe K A_1 is now p-iven hv 


, K+ i, with z K = 


X k,~ 1-i 


•” ' 2 -“2 

The Wasserstein-like metric (and its corresponding norm) needs to be adapted as well: the scalar 
product (•, -y s is replaced by 

<v,w> ? = Yi S k v kW k and ||v|| 5 = <v, v)g. 
fee I° K 

Hence the metric gradient V^/(x) e R A+1 of a function / : —> R at x e is given by 

[V,T(x)]fe = ^*/(x). 

Otherwise, we proceed as before: the entropy is discretized by restriction, and the discretized 
information functional is the self-dissipation of the discretized entropy. Explicitly, the resulting 
fully discrete gradient flow equation 


'-A 


= -VrF a , A (x£) 


attains the form 


X u Xi 


2a 


(2a + 1) 2 4 


[(^ +| )“ + i[D|(z”)“ + l] fe+| ,-(^_i)“ + i[D|(z”)“H] fe _r] +Ax^, 

(67) 


with [ D |z“ + 5 ]fc _ i : = (z::? - 2 a+cm 

4.2. Implementation. To guarantuees the existence of an initial vector e $s, which ’’reaches” 
any mass point of u°, i.e. [.Tq, x a K \ c supp(n°), one has to consider initial density functions u° 
with compact support. 

Starting from the initial condition x^, the fully discrete solution is calculated inductively by 
solving the implicit Euler scheme (671 for x^, given x^ -1 . In each time step, a damped Newton 
iteration is performed, with the solution from the previous time step as initial guess. 

4.3. Experiment I — Exponential decay rates. As a first numerical experiment, we want to 
analyse the rate of decay in case of positive confinement A = 5, using a = 1. For that purpose, 
consider the initial density function 


u° = 


0.25| sin(a;)| • (0.5 + ILexj^)), x e [— n, 7r], 
0, esle 


( 68 ) 


Figure [3] shows the evolution of the discrete density ua at times t = 0.05,0.1,0.15,0.175,0.225, 
using K = 200. The two initially seperated clusters quickly merge, and finally changes the shape 
towards a Barenblatt-profilc (doted line). 

The exponential decay of the entropies H aiA and F ajA along the solution can be seen in figure 
|T]/left for K = 25, 50,100, 200, where we observed the evolution for t e [0,0.8]. Note that we write 
H a ,\(t) = H ajA (x^) and F a< \(t) = F a>A (x^) for t e (t n _ and set ifo, a( 0) = H a>A (x^) and 
F fliA (0) = F ajA (x^). As the picture shows, the rate of decay does not really depent on the choice 
of AT, in fact the curves lie de facto on the top of each other. Furthermore, the curves are bounded 
from above by 0H a>A (u°) — H a A (b ttiA )) exp(—2A£) and (F a , )A (u°) — ^^(ba^)) exp(— 2\t) at any 
time, respectively, as (16) & (17) from Theorem [2] postulate. One can even recognize, that the 
decay rates are even bigger at the beginning, until the moment t = ... when ua finishes its 









LONG-TIME BEHAVIOUR FOR A FAMILY OF FOURTH ORDER 


23 



Time t 



Figure 1. Left: Numerically observed decay of H ai \(t) — and F a ,\(t) — 
F™ along a time period of t e [0, 0.8], using K = 25,50,100,200, in compar¬ 
ison to the upper bounds ("H Qj \{u°) — exp(—2 Xt) and (Fan'll 0 ) — 

*T*a,A(bo;,A)) exp(—2At), respectively. Right: Convergence of discrete minimizers 
w™ ln with a rate of A' -15 . 


’’fusion” to one single Barenblatt-like curve. After that, the solution’s evolution mainly consists 
of a transveral shift towards the stationary solution b aj ,\> which is reflected by a henceforth 
constant rate of approximately — 2A. 

Moreover, figure fright pictures the converence of rt“ ln towards b Qj ,\. We used several val¬ 
ues for the spatial discretization, K = 25,50,100,200,400,800, and plotted the L 2 -error. The 
observed rate of convergence is AT -15 . 


4.4. Experiment II — Self-similar solutions. A very interesting consequence of section [312 
is, that the existence of self-similiar solutions bequeath from the continuous to the discrete case. 
In more detail, this means the following: Set A = 0 and define for t e [0, +oo) 


b a,o(T •) := 3fi(t)b Q ,i, with R(t) := (l + (2a + 3)t) 


l/(2a+3) 


(69) 


then b* 0 is a solution of the continuous problem (551 with u° = b* 
solutions to 


Q,0 


(0, •). In the discrete setting, 
with A = 0 are inductively given by an initial vector b^ a 0 with corresponding 

density = uy[b^ 0 ] that approaches b* o (0, •), and t>A,a,o = ^s? b A,aO wb; h S'" defined as 
in (631, for further n = 1,2 

As figure [2] shows, the resulting sequence of densities ua (black lines) approaches the con¬ 
tinuous solution b* 0 of (69) (red lines) astonishingly well, even if the discretization parameters 
are choosen quite rough. In this specific case we used K = 50 and r = 10” 3 . The discrete and 
continuous solutions are evaluated at times t = 0,0.1,1,10,100 
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Space x 


Figure 2. Snapshots of the densities b* 0 (i, •) (red lines) and «a (black lines) 
for the initial condition b* o (0, •) at times t = 0 and t = 0.1 • 10®, i = 0,..., 3, 
using I\ = 50 grid points and the time step size r = 10 -3 . 



-4 -2 0 2 4 -4 -2 0 2 4 -4 -2 0 2 4 

Space x Space x Space x 



-4 -2 0 2 4 -4 -2 0 2 4 -4 -2 0 2 4 

Space x Space x Space x 


Figure 3. Evolution of a discrete solution u&, evaluated at different times 
t = 0,0.05,0.1,0.15,0.175,0.25 (from top left to bottom right). The red line is 
the corresponding Barenblatt-profile b aj A- 
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